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EMTEX 



1 . The EMTEX beamline is part of the Universal Linear Accelerator (UNILAC) 
at GSI. 

2. EMTEX stands for EMittance Transfer Experiment. The aim is to 
demonstrate emittance transfer, e.g. from vertical to horizontal. 

3. Positive deuterium ions are injected into EMTEX. A carbon foil in the 
beamline strips electrons from the ions. 

4. Emittance transfer is achieved by placing the foil at the centre of a solenoid, 
together with a suitable design of the beamline magnets. 
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A requirement 



To measure: 



1 . The four-dimensional symmetric beam matrix 

" (xx) (xx f ) (vy) (*y') " 
(x f x) (x f x f ) (x f y) (x f y f ) 
(yx) (yx f ) {yy) (y\ J ) 

_<A> Wm <y f y) <y'y f )l 



c = 



2. The 4D phase space distribution. 



Setup: 



To use a vertically focussing quadrupole, a 
horizontally focussing quadrupole and a 
scintillating screen. 




EMTEX beamline 



For this feasibility study, we assume that the 
is a screen halfway between the normal and 
skew triplets. 

We assume that the 4D measurement will 
make use of this screen and the two 
quadrupoles just before it. 



Vertically focussing 

Horizontally focussing- 
Assume a screen here" 



C. Xiao, et al (2013) 
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Foil stripping 



At the foil, charges are stripped from deuterium ions with energy 11 .4 
MeV/u. 

In publications, different notations have been used for these ions: 
At injection: 



After stripping: 




where D is a deuterium atom with 1 proton and 1 neutron. 
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Problems and suggestions 



Problems 

1 . There is a standard method for measuring 2D beam distribution matrix, but it is not 
clear how to measure a 4D beam distribution matrix. 

2. Recent simulation work has shown that 4D phase space measurement is possible, 
but this has not yet been demonstrated by experiments. 



Suggestions 

1 . Standard 2D beam distribution matrix method may be generalised to 4D. There 
many unknown parameters to find. We need to check if it is robust. 

2. 4D phase space measurement method could be simulated for EMTEX to 
determine its feasibility. (The 4D beam distribution matrix may also be determined 
from the phase space distribution.) 
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2D beam distribution matrix measurement 



One way is to use a quadrupole scan setup: 

1 . Take a photograph of the beam image at B. 

2. Change quadrupole strength. 

3. Take another photo. 

4. Repeat steps 2 and 3 for 1 00 - 200 times. 

5. Then find the beam distribution matrix from the beam images. 

Find beam distribution matrix here 



Electron bunch 





y 



^ screen 
beam image 



Finding the 2D beam distribution matrix 



The beam distribution matrix at the entrance to the quadrupole is 



C = 



r , 2 
x 



A A 



V 



JCJC 



JCJC 

x 



J 



The i th transfer matrix from entrance of quadrupole to screen is 



AT, = 



(M n . M 



\ 



12, i 



yM 2li M 22i j 



The beam distribution matrix at the screen is 



C = 



// 2 
X 



A \ 



v 



X 



'2 



ij 



They are related by 



C = MCM) 



Finding the 2D beam distribution matrix 



Expand and pick out the <x 2 > term: 



* 2 ) . = M i u (x 2 ) + 2M l U M U . (xx) + M 12 , 2 (x 



'2 



If n images are measured, then there are n such equations. Matrix M can be 
computed from the quadrupole parameters. 

We can then find <x 2 >, <xx'> and <x' 2 > that best fit the equations. 



beam image 



C\J 



E 



CM 
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beam distribution matrix 



C = 



r, 2 
x 



V 



XX 

'2 

X 



J 



X 



k^fm' 1 ) 
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But how to find 4D beam distribution matrix? 



The beam distribution matrix at the entrance to the quadrupole is 



c = 



ft 2 
X 



XX 



v 



xy 

xy' 



XX 

'2 
X 

xy 
xy 



xy 
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yy' 



A A 

xy 

xY) 
yy 

'2 
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j 



The i th transfer matrix from entrance of quadrupole to screen is 



m = 



r M xu 


M 124 


0 


0 ^ 




M m 


0 
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The beam distribution matrix at the screen is 

n 2 
x 



C = 



v 



xy 
xy 



XX 

1 

X 

xy 
xj 



xy 
xj 



xy 



A A 



yy 



y 



'2 



ij 



They are again related by 



G = M CM 



ii 



One possible way 



The quantities that can be obtained from the beam images are <x 2 > h <y 2 >j and <xy>j. 
Try writing down the expanded equations: 



(xy). = M lu M 33J (xy) + M n .M 33 .{xy) + M nj M 34J (xy) + M l2 .M 34J (x'y') 



The right hand sides contain all 10 independent elements in the 4D beam distribution 



The 1 st and 3 rd equations can be solved separately using the 2D method for <x 2 >, 
<xx'>, <x' 2 >, <y 2 >, <yy'> and <y' 2 >. 

The 2 nd equation contains 4 unknowns <xy>, <x'y>, <xy'> and <x'y'>- 

In principal, they can be solved by finding a set of values that best fits the above 
equations, but we should do careful simulation to check if this method is reliable. 12 





matrix. 



Phase Space Measurement 



1 . The technique described in the previous slide gives the second moments. It is 
also interesting to look at the phase space. 

2. There is a standard method using tomography to measure 2D phase space 
distribution at a point in a beamline. 

3. This method requires a screen to record beam profile, some way to vary the 
transfer matrix before the screen, and a computer code to process the data. 

4. The standard method can only work for 2D. Recently, we have developed a 
method to measure 4D phase space distribution. 

5. We have demonstrated using simulation the 4D distribution is possible in the 
ALICE accelerator at Daresbury using a vertically and a horizontally focussing 
quadrupole, and a screen. 

6. However, the method has not yet been tested in experiments. We shall show 
by simulation that it is also feasible in EMTEX. 
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2D phase space measurement ? 



In the hospital, tomography is used to image the internal organs. It works by taking 
X ray projections of the patient from different angles. 




http://www.fda.gov/radiation- 

emittingproducts/radiationemittingproductsandprocedures/medicalimaging/medicalx- 
rays/ucmll5318.htm 

If we can take projections of a 2D phase space from different angles, then we can 
also do tomography on the phase space. 

14 



Phase space tomography 



Experimental setup 

1 . We want to find the 2D phase space at a point A in a beamline. 

2. The phase space at a screen B can be rotated if we change the transfer matrix 
between A and B (e.g. by adjusting a quadrupole in between). 

3. We can then use a standard method to find the phase space distribution. 



A 



y 




Screen 



x 



A 




B 



quadrupole 
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Phase space tomography 



Theory: 1 

2 
3 
4 
5 



Projection of screen image in x direction = projection of the phase space. 
This is related to projection in rotated direction at A (reconstruction location) 
The rotation angle can be obtained from the transfer matrix. 
By changing the transfer matrix, we can get projections for different angles. 
Then reconstruct using standard tomographic technique. 

Screen 




A 



B 



x 



x 




X 



B 



\ X B J 



M u M 



12 



A 



\M 2X M 22 J 



\ X A J 



x 




Filtered Back Projection (FBP) 




To reconstruct original distribution from projections: 



1 

2 
3 



Fourier transform each projection, filter, invert transform 
Spread over the phase space in perpendicular direction 
Repeat for each angle and add over previous one. 



FFT, filter, invert 



spread 



1 projection 



4 projections 



next 




Kak and Slaney (1999) 



simulations 



iections 




(close to original) 
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Example: Measurement at Daresbury 



ALICE tomography 
section 



1 . Take a photograph of the beam image at YAG-02. 

2. Change quadrupole strength of QUAD-07. 

3. Take another photo. 

4. Repeat steps 2 and 3 for 1 00 - 200 times. 

5. Use Filtered Back Projection (FBP) to reconstruct distribution at A. 
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4D phase space measurement 



Standard phase space tomography works for 2D phase space (x,x') 5 (y,y') or (z,5) only. 



Is it possible to measure and reconstruct a full 4D phase space like (x,x', y,y')? 
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Simulations 



19 



4D phase space measurement 



Reconstruction 
location 



screen 



D 

QUAID-06 
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ALICE tomography 
section 



A possible setup: 

1 . Suppose that a horizontally focussing quadrupole (F) can rotate 
in (x,x') and a vertically focussing one (D) in (y,y'). 

2. For each (y,y') angle, record the images for different (x,x') 
angles. 

3. Repeat this for different (y,y') angles. 

4. The full set of recorded images can then be processed using a 
series of 2D reconstructions to reconstruct the 4D phase space. 



Details, equations and algorithm are given in: 



K. M. Hock and A. Wolski. "Tomographic Reconstruction of the Full 
4D Transverse Phase Space", submitted to NIMA (2013). 
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Simulation of a measurement 



Using 2 quadrupoles in ALICE, range of angles limited to 150° in both x and y. 



Reconstructed 4D distribution still shows good agreement. 
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4D phase space measurement 



Standard phase space tomography works for 2D phase space (x,x') 5 (y,y') or (z,5) only. 



Is it possible to measure and reconstruct a full 4D phase space like (x,x', y,y')? 



Projections from 
a hypothetical 4D 
distribution 
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Simulations 
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4D theory: From 4D to 2D 



Define 4D phase space distribution at A (reconstruction location) as number density of 
electrons per unit phase space volume: 

Suppose that this is mapped to B, the screen, as 

f B (x,x\y, y') 

Screen image is projection of 4D phase space to 2D 
(x,y) plane of screen: 



l(x, y) — 




f B {x y x\y 7 y')dx'dy' 



Think of this in 2 steps: a projection to a 3D (x,x',y) phase space 

Q y {x,x') = J f(x,x\ y, y f )dy'. 

followed by a projection of this 3D phase space to the 2D screen 

I(x,y) = jg y (x,x')dx' 




Visualise the 3D phase space with an x' axis perpendicular to the screen. 
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4D theory: From 2D to 3D 



Then we can think of a horizontal line on the screen as a projection from the 3D phase 
space in the x direction - for a given y. 



If we can somehow rotate the (x,x') space about the y 
axis, we can collect a set of these projections and 
reconstruct the distribution for each y. Then if we repeat 
this for every y, we would have reconstructed the (x,x',y) 
phase space. 




From 2D phase space reconstruction, we already know 
that we can rotate the (x,x') space using a horizontally 
focussing quadrupole. So we know how to reconstruct 
(x,x',y). 

Recall that this 3D phase space is a projection from the 4D phase space: 

f(x,x',y, y f )dy l . 




Think of a vertical line in this 3D phase space for a given (x,x'). This is actually a 
projection from the 4D phase space in the y direction. 



4D theory: From 3D to 4D 



Since the vertical line through (x,x') is a projection in y direction, we may be able to 
reconstruct (y,y') for each (x,x') - if we could rotate the (y,y') phase space. 



Practical issues 

The 4D simulation results use 1 horizontally focussing and 1 vertically focussing 
quadrupole. These give limited ranges of angles of 150° in x and y, resulting in slight 
distortion and higher background. 

Use of more quadrupoles can potentially increase the range of angles to full 180°. 



Since we can reconstruct (y,y') for each (x,x') , the full 
4D phase space can now be reconstructed. 



From 2D phase space reconstruction, we know that we 
can - by using vertically focussing quadrupole. 




(x,x') 



x 



X 
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Feasibility of 4D tomography at EMTEX 



1 . The previous slides show the feasibility of 4D tomography in an electron beam 
accelerator at Daresbury. 

2. Can the optics in the beamline be adjusted to provide a sufficient range of 
rotation angle? 

3. Another questions is the beam size, This expands during vertical focussing. Is 
the screen big enough? 

4. To answer these questions, we simulate the motion of 30000 D 3 + down the 
EMTEX beamline. 
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EMTEX beamline 



For this feasibility study, we assume that the 
is a screen halfway between the normal and 
skew triplets. 

We assume that the 4D measurement will 
make use of this screen and the two 
quadrupoles just before it. 



Vertically focussing 

Horizontally focussing 
Assume a screen here 



C. Xiao, et al (2013) 
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Simulation method 



1 . Coordinates (x,x',y,y') of a random, 4D Gaussian distribution of 30000 
particles is generated, with beam parameters given on the previous slide. 

2. Each particle has a mass of 6m u , energy of 6x11 .4 MeV and charge of +e. 

3. It is propagated from start to end of the lattice described on the previous 
slide. 

4. Linear beam dynamic is assumed. So 4x4 transfer matrices are used for 
each element of the lattice, except for the foil. 

5. At the foil: 



(b) There is an increase in angular spread of A§ = 0.226 mrad in both x 
and y directions. So the standard deviation of x' and y' in the 4D 
distribution are both increased by A<|>; 




becomes 3 particle; ® + ® + ® + 



(c) Each particle now has a mass of 2m u , energy of 2x11 .4 MeV and 
charge of +e. 
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Simulation of 



10 cm x 10 cm 
screens assumec 



start 



halfway between 
normal and 
skew triplets 



end 



1 
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profiles 
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Check simulation against publication 



C. Xiao, et al (2013) 

TABLE IL The beam parameters at the entrance and exit of the 
EMTEX beam line. An initial w ate r- bag distribution was as- 
sumed for ihe tracking simulations. 



Parameters 
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SIMULATION. Parameters computed from second moments of 
beam distribution. 



Parameters 


Entrance 


Exit 


Px/P y [m] 

8 x /e y [mm. mrad] 


-1 .23/-2.28 
22.10/15.02 
0.509/0.509 


0.19/0.04 
6.64/6.46 
0.306/1.163 



Some differences 
Could be due to 
foil simulation. 
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Projection angles at entrance face of D quad just before screen 

D F screen 




Scan F or D quad just before screen, keeping the strength of the other quad fixed 



Scan D quad 
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Scan F quad 
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Rotation Angles 




1 . Rotation angle is determined by transfer matrix. 



2. For good reconstruction of phase space, it should cover the full range of 
180°. 



3. For 4D measurement, the full vertical and horizontal range should ideally be 
accessible for any combination of vertical or horizontal angle. 

4. Results in previous slide show that the accessible ranges of angles are 
reasonably close to 180° for most but not all angle combinations. 



Accessible ranges of angles 



1 . The vertical 0 V and horizontal 6 h phase space rotation angles for possible 
combinations of F and D quad strengths are calculated and plotted. 

2. The resulting graph shows that at least 150° x 150° is accessible. An earlier 
slide has shown that this is enough for a 4D measurement. 



180 



This region is not accessible. 
It is a limitation of the two 
quadrupole setup. Using more 
quadrupoles may help. 




e h (deg) 



1 50° x 1 50° 
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Largest beam spots in y and x directions 

D F screen 



10 cm x 10 cm 
screen assumed 




D quad 0 T/m, F quad 8.421 T/m 



D quad -9.431 T/m, F quad 0 T/m 





x (cm) 



x (cm) 
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4D tomography at EMTEX 



1 . Simulation shows that with the bipolar quadrupoles, the ranges of vertical / 
horizontal angles cover at least a 1 50° x 1 50°, which is sufficient. 

2. It is reasonable to think that it is also sufficient for beam distribution matrix 
measurement since tomography requires more data, but this needs further 
study. 

3. Typical screen size is large enough to capture the defocussed beam images 
for the assumed starting beam. 

4. The above results suggest that 4D phase space measurement at EMTEX is 
feasible. 
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Conclusion 



1 . The requirement to carry out measurements of 4D beam distribution matrix 
and 4D phase space at EMTEX is described. 

2. A standard method for 2D beam distribution matrix measurement is 
reviewed. 

3. A4D beam distribution matrix measurement is proposed by generalising the 
2D method. 

4. 2D phase space measurement using tomography is reviewed. 

5. Recent development of 4D phase space measurement is reviewed. 

6. 4D measurement at EMTEX is studied by simulation and shown to be 
feasible. 



References 



C. Xiao, et al, "Single-knob beam line for transverse emittance partitioning", Physical 
Review Special Topics - Accelerators and Beams, 16, 044201 (2013). 

H. Wiedemann, Particle Accelerator Physics, 3 rd ed., Springer- Verlag, Berlin (2007), pp. 
164-166. 

M. G. Ibison, et al, "ALICE Tomography Section: Measurements and Analysis," Journal 
of Instrumentation, 7, P04016 (2012). 

K. M. Hock and A. Wolski, "Tomographic Reconstruction of the Full 4D Transverse 
Phase Space", submitted to Nuclear Instruments and Methods A (2013). 



37 



Thank you 
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